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Abstract. We present a detailed symmetry analysis of the degeneracy lifting due to 
higher order spin exchanges in the pyrochlore lattice in applied magnetic field. Under 
the assumption of the 4 sublattice ordering, the criteria for a stable half- magnetization 
plateau are deduced. The higher order exchange terms may originate from spin-lattice 
coupling, or can describe quantum and thermal fluctuations. 
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1. Introduction 

The very special feature of the antiferromagnetic nearest-neighbour Heisenberg model 
on a pyrochlore lattice is that it does not order magnetically at any temperature [UI21- 
This is a consequence of the high (continuous) degeneracy of its ground state manifold, 
which for 0(3) spins is too great even for order-from-disorder effects to be effective. 
None the less, the vast majority of magnetic materials with a pyrochlore lattice do 
order magnetically. In these compounds, magnetic order is typically accompanied by 
a structural transition which lowers the crystal symmetry and so lifts the ground state 
degeneracy. This "order by distortion" mechanism [31 E] is a natural consequence of 
the magnetoelastic interactions present in almost all magnets. Examples may include 
YMn 2 E], Y 2 M 2 7 but to date its most striking realization is in the high field 
properties of the chromium spinels CdCr 2 04 |Zj and HgCr 2 04 [Hj. 

In these compounds a dramatic half-magnetization plateau is accompanied by a 
collossal magnetostriction and change of crystal symmetry. The basic physics of this 
magnetization plateau can be understood by extending the simplest four-sublattice 
"order by distortion" model of the pyrochlore lattice [3 E] to finite magnetic field [Oj. 
Here we develop a complete symmetry analysis of the most general four-sublattice order, 
and give additional details of the new phases which result when the assumptions made 
in |Hj are relaxed to allow for more general spin-lattice interactions. 

2. The effective model 

We start by assuming a simplest microscopic form of magnetoelastic coupling, in which 
the strength of exchange interactions between a pair of nearest neighbour classical spins 
depends linearly on the change in length 5ij of the associated bond 



Here a is the spin-lattice coupling, and K an elastic coupling constant. We choose 
the z-axis to coincide with the direction of applied magnetic field. As the experiments 
point to an ordered state, we further assume that the simplest, four-sublattice order 
is stabilized, e.g. by longer range exchange interactions. Since we are interested in 
static quantities, we can then safely reduce the problem to that of a single tetrahedron 
embedded in a pyrochlore lattice \. 

The symmetry group of this embedded tetrahedron is 7^, which has 24 elements 
and 5 irreps. Bond variables such as 5y and Sj • Sj transform according to the Ai, E 

| In general, the changes in bond length 5i_j are not independent variables, but rather implicit functions 
of the positions of all the magnetic ions. However, since in this special four-sublattice case the lattice 
is subject to a uniform distortion (an affine transformation), the lengths of the 6 bonds making up the 
tetrahedron can be treated as independent quantities. 
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In this approach we assume that the the spin-space and the real-space are decoupled, 
i.e. that the L-S coupling and crystalline anisotropics are negligible. 




A T 2 ,1 A T 2 ,2 A T 2 ,3 h/J 



Figure 1. (a) Schematic representation of the different irreps of the tetrahedron. 
Hollow red lines have positive weight; green lines have negative weight; thin black lines 
have zero weight, (b) Magnetization curve for 6ai = 0.1, &e = 0.1 and bj 2 — 0.2. The 
spin configurations, together with the irrep they belong to are indicated. 



Under these assumptions, the Hamiltonian of the embedded tetrahedron is given 

by 

H = 2v / 6JAa 1 - 4/iM - 2J (a Al A Al p Al + a E A E • p E + a-r 2 A T2 • p T J 

+ (K Al pl i+ K E pl + K T2 p 2 T2 ) . (3) 

M is the average magnetization per site. The energy minima are found for pr = 
(aRj/KR)AR, so the energy above becomes 

H = 2J (VQA Al - b Al A 2 Al - b E A 2 E - b T2 A 2 T ^ - AhM (4) 

with 6 R = JolrI '(2Kr). Actually, for the simple model defined by Equation (JTJ), all the 
couplings turn out to be equal: b Al = 6 E = bj 2 , and this is the case we have considered 
in Ref. This implied that only biquadratic terms of the form (Sj ■ Sj) 2 were present 
in the Hamiltonian. Once we allow nonequal couplings, we get three and four-site terms 
in the effective Hamiltonian of the form 

H = 2 J ^ • - h St - J -(b Al + 2b E + 36 T2 ) $^(S 4 • S,) 2 

i,3 i i,j 

2 T 

- T (b Al -b E )J2(Si-Si)(Si-S k ) 
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- ^(6 Al + 26 E - Sb T2 ) (Si ■ Sj)(Sk ■ S,) , (5) 

i,j,k,l 

where the summations are carried out such that every possible combination of spins in 
the tetrahedron appears exactly once, furthermore i 7^ j 7^ k 7^ I. 

This result is of a very general nature: Equation (jlj) is the most general Hamiltonian 
at fourth order in spin operators describing a 4 sublattice long-range ordered state on 
the pyrochlore or fcc-lattice. As such, it is independent of the actual origin of the 
coupling, which may equally be due to quantum and thermal fluctuations, or higher 
order exchange processes JU] ■ It could have been derived on purely phenomenologically 
grounds, as the sum of the group invariants — A Al = 8(M 2 — l/A)/y/E (and its powers) 
is itself invariant while we can construct the quadratic forms = A| 1 + A| 2 and 
Aj 2 = Aj 2 1 + Aj 2 2 + Aj 2 3 from the E and T 2 irreps. There are six cubic invariants; 
three trivial ones A A , A Al A|, and A Al Aj 2 , and three nontrivial ones A| 1 — 3Ae,iA| 2 , 
At 2 ,iAt 2 ,2At 2 ,3 and 

/ 1 . Fx \ I 1 . Fx \ ^ 




Ae,iA^, 3 • (6) 

These can also, in principle, be added to the Hamiltonian, but resulting proliferation of 
parameters would make the problem intractable, and for that reason we neglect them 
here. All of the invarients considered above are constructed from two-spin dot products: 
at 6th order more complicated spin expression may appear. 



3. Stability analysis of the collinear states 

Both the plateau and the fully polarized states are realized with collinear spin 
configurations, with spins pointing along the magnetic field. At certain critical values 
of magnetic field, these configurations cease to be favorable, and the spins cant to 
optimize their energy. This can be a continuous transition, when the magnetization 
smoothly changes from the fractional value, or a first order transition, in which case the 
magnetization jumps. Here we introduce a simple effective Hamiltonian to treat this 
effect, first for the fully polarized state, and then for the plateau. 



3.1. Instability of the fully polarized state 

Small fluctuations about the fully polarized state can be characterized by 

S, ('(rCry 1 '/J ^j) • (7) 

where rjj and £j are small. The rjj and £j are site variables, which transform according 
to the Ai and T2 irreps§. Furthermore, since the model is isotropic, the magnetization 
always points along the magnetic field, J2j Vj = J2j £j = 0- This means that the A 1 

§ The bond variables, which are "squares" of site variables, can transform as Ai, E and T2 (as the 
product J 2 ® T 2 = Ai ® E © Ti © T 2 ). 
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irrep is zero, and we need to worry about the T2 irrep only. Apart from the space- group 
symmetry, the magnetic field reduces the 0(3) spin-symmetry to U(l), the rotation 
of spins around the axis parallel to the magnetic field. Instead of the real rjj and £j 
quantities it is useful to introduce the complex rjj + i£j quantity: 

Vi + ifi =^1+^2 + ^3, V3 + = -Hi + Hi ~ H3 , ^ 
V2 + i£,i = Hi - Hi - H3 , V4 + i& = ~Hi ~ Hi + Hz > 

where the complex /x = [fi\, Hz) transforms according to the 3-dimensional T2 irrep. 

Assuming a second order transition, the fully polarized state becomes unstable at 
a critical field of 

hp = 8J(1 - 26 Al ) ■ (9) 

Expanding to fourth order in fi at this critical field, the energy per spin of the resulting 
canted state reads : 

E = E F + hji - M|/x| 2 + ^ (3 - 226 Al - 8b E ) |/x| 4 

+ 8J(b E - b T2 ) (\hi\ 2 \H2\ 2 + |/ii| 2 |/i3| 2 + \H2\ 2 \H3\ 2 ) 

1 71 / 2 — 2 1 — 2 2 1 2 — 2 ,—22, 2—2 1 —2 2\ /inN 

-4J6j 2 {H1H2 + H1H2 + H1H3 + H1H3 + H2H3 + H2H3) , (10) 

where Ep = 3J(1 — b^) — h is the energy of a spin aligned with the magnetic field. 
Minimizing this energy as a function of the coupling parameters we obtain the results 
shown in figure Efa). In clockwise direction, we observe four phases: (i) The coplanar 
2:2 canted phase of E symmetry with fi cx (1, 0, 0) (and permutations) when b E > and 
&e > 26j 2 ; (h) The T2 symmetric coplanar 3:1 canted state with = | /X2 1 = | /X3 1 and 
equal arguments modulo 7r for bj 2 > and b E < 2bj 2 ; (iii) T2 symmetric umbrella like 
configuration with fi cx (1, e"^, e~^) when bj 2 < and b E < 2bj 2 ; (iv) An E symmetric 
1:1:1:1 umbrella like configuration with fi oc (1, ±2, 0) when b E < and b E > 2bj 2 . 

The transition changes from a continuous into a first order one when the prefactor 
of the quartic term becomes negative. In the 2:2 phase, this happens for 8&e > 3 — 226 Al , 
while in the 3:1 state for 16bj 2 > 3 — 22bp Kl . Similarly first order transition takes place 
when 22&Ai > 3 — 4bj 2 in case (iii) and 22b ^ > 3 — 2b E in case (iv). 



3.2. Instability of the plateau 

The analysis can be repeated for the stability of the plateau state: we assume that the 
spin at site 1 is down, Si = ^771, £1, — \/l — n\ — , while the others (sites 2, 3 and 
4) are up. Instabilities occur at 

hi = 4J{l-4b T2 ) and h u = 4J(1 + 26j 2 ) , (11) 

i.e. the plateau is stable for hi < h < h u , assuming a continuous phase transition, and 
its width is 24bj 2 . Clearly first order transitions, if they occur, will shrink the magnitude 
of the plateau. Let us consider what happens at the upper edge of the plateau. The 
selection of the position of the down pointing spin breaks the symmetry in the [i vector, 
and in our case, the soft mode is the /x oc (1, 1, 1) direction. It is convenient to introduce 




Figure 2. Instability of the (a) fully polarized state, and of the plateau at the (b) 
upper critical filed h u and (c) lower critical field hi . The magnetic field is perpendicular 
to the plane, and the red arrows shows the spin canting (f]j,£j) of the j-th spin. Green 
dots and dashed arrows represent the down pointing spin in (b) and (c), while the 
shaded grey region where the plateau is not present. 



the /J = (/ii + /i 2 + [i3)/\/3, and the fj,f = (2/ii - /i 2 - [13) and ^ = (/i 2 — 
orthogonal to the fjfi, the energy up to fourth order in fj, (keeping the relevant quartic 
term only) reads 

E = E P -~(h- h u )\^\ 2 + 12Jb T2 (|^| 2 + |/4| 2 ) + ^(3 - 46 Al + 26 T2 )|/i ll | 4 • (12) 

Here E P = —h/2 — 3Jbj 2 is the energy of the plateau. The instability is toward the 
3:1 coplanar state, with the up spins having all equal canting angles. The continuous 
transition becomes first order once 3 — 4b ^ + 2bj 2 < 0. This is the same 3:1 canted 
coplanar state found as the instability of the fully polarized solution, and there is no 
reason to suppose that the two states are not connected with each other - in other words, 
there is no symmetry requirement for an additional phase transition when moving from 
the plateau to the fully polarized state ||. 

Next, let us consider what happens at the lower edge of the plateau: the energy is 

E = E P - - h) (l/^l 2 + I/4I 2 ) + 6&t>"| 2 

2J 2 
+ y (4& E + bj 2 ) \n x M 2 - M X M2 I 

+ f (3 - 46 Al - 8b E - 246 T2 ) (|^| 2 + |/4| 2 ) 2 (13) 

and, contrast to what happened at the upper critical field h u , the softening occurs in the 
two modes perpendicular to /J. The l/x^/i^ — fii^l term selects the relative phase of 
the fii and /i^: /ii" = ±^2") an d a T 2 symmetric umbrella-like canting is realized for 
|| This need not be the case if we consider ordered states with more than four spins in the unit cell. 
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46 E + bj 2 < 0, while for 46 E + bj 2 > the instability is against a degenerate coplanar 
canting (only the phase of the /i^ and is fixed; their relative amplitude ^2/^1 * s 
not). In fact we need to go to 6th order to select the 2:1:1 state of mixed E and T 2 
symmetry previously found numerically j^j. This has three solutions: (/ii, /i 2 ) (1,0), 
or (—1/2, ±\/3/2). These transitions become first order when 226j 2 > 3 — 4&Ai for the 
umbrella state, and when 86 E + 246j 2 > 3 — 46m f° r 2:1:1 state (see figure I2tc))- 



4. Extrema of the invariants and magnetization curves 

Any given magnetization can be realized by a large number of (continuously degenerate) 
spin configurations. Within our model, the energy of these states is uniquely determined 
by the second order invarients A R . Therefore, by mapping out the extrema of these 
invariants, we exhaust all possible ground states selected by four spin interactions 6r. 
Doing so enables us not only to confirm the results of the stability analysis above, but 
to make some further observations about the magnetization process M(h) and possible 
first order transitions. 

Let us first consider the limiting case Ay 2 = 0. Then, from (J2J), the exchange 
energy on the opposing bonds of the tetrahedron is equal. This condition is satisfied 
by spins with equal S z components, and the (S x ,S y ) = (±1,0) and (±cos0, ±sin0). 
The maximum of value of A| is obtained for a 2:2 canted state with = 0, while the 
minimum is obtained for a fully symmetric 1:1:1:1 canting with = tt/2 (see figure EJ), 
such that 

l(l-M 2 ) 2 <A 2 <^(l-M 2 ) 2 . (14) 

Next, let us consider the case of A E = 0. The maximum value for Ay 2 = 6 is realized 
when all the Aj 2i j take their maximal value Aj 2i! = V2 in the collinear uuud (three up-, 
one down-spin) configuration with M = 1/2. For M > 1/2, Aj 2 is maximal for the 
coplanar 3:1 canted state. For M < 1/2 the maximal state is a a three-dimensional 
inverted umbrella type configuration (Sj = (sin?? cos[j27r/3], suit? sin[j27r/3], cosfi) for 
j = 1, 2, 3 and S4 = (0, 0, —1), while the minimal is obtained by inverting S4 to obtain 
a folded umbrella configuration, and 

^M 2 fl-M) 2 <A 2 <{ K 1 -^) 2 > ifV2<M<l; 
3 M{1 M) - Aj2 -\ fM 2 (l + M) 2 , ifO<M<l/2. (15) 

However, a slightly better energy can be obtained by allowing finite values for A 2 :, in 
which case the extremal point corresponds to coplanar 2:1:1 canted state with a relatively 
complex analytic form. 

The existence and stability of the plateau is related to the cusp of Aj 2 at M = 1/2. 
As a singular point, it will minimize the energy for a finite window of the magnetic field 
when 2bj 2 > max(6E,0), also allowing for negative values of 6 E - This is in agreement 
with the results of the stability analysis presented above. The extrema in the (A E , Aj 2 ) 
plane can clearly be related to the nature of instability of the plateau and the fully 
polarized state. 
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Figure 3. (a) The extremal values of A| and Aj 2 as a function of the magnetization 
M. (b) The cut at M — 3/8. The degenerate manifold maps into the shaded region, 
and the border is outlined in blue for coplanar and red for non-coplanar states. The 
spin configurations are shown from the direction of the magnetic field, with red arrows 
canted in the direction of the field, and the green dots represent spins pointing opposite 
to the field. 



However, we can also see that for > max(2&r 2 > 0) the T 2 order is not realized at all 
(here we allow for negative values of bj 2 ). The system remains in the coplanar 2:2 phase 
for all values of applied field up to saturation. The half magnetization plateau is also 
absent when both and bj 2 are negative: the ground state in this case is a noncoplanar 
spin configuration. It is also possible to achieve a T 2 type solution for small field, and 
an E type configuration for high field, with a first order transition between them. 

In figure IHb) we show the magnetization curve for a physically relevant set of 
parameters. To facilitate comparison with the experimentally measured curves, it is 
useful to extract the characteristic features of the magnetization curve. The energy per 
site of the 2:2 canted E symmetry state is: 

E 2:2 = -KM - J(l - AM 2 ) - il - AM 2 ) 2 b Al - ^(1 - M 2 ) 2 b E (16) 
with zero field susceptibility 

SM 3 1 (17) 



6h 



h=0 8J3 + 26 Al +46 E 



Similarly, the energy of the 3:1 canted state for M > 1/2 is : 

E 3a = -hM - J(l - AM 2 ) - ^(1 - AM 2 ) 2 b Al - ^(1 - M 2 ) 2 bj 2 , (18) 

with differential susceptibility at the upper critical field: 
5M 3 1 



Sh 



= - fig) 

h =4j+s bT2 8J3-46 Al + 26 T2 ' 

If we continue the tangent of the magnetization curve at the upper critical field h u of 
the plateau, it will cut the M = axis at h = 16J(b Al + bj 2 )/3 (Ahi in figure HJb)). 
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Similarly, if we continue the linear magnetization at low field, it will intersect the plateau 
at a distance Ah 2 = 8J(3bj 2 — — 26e)/3 from h u . These values are easily extracted 
from the experimentally measured magnetization curves, and together with the upper 
critical field for the plateau, ([lip, allow to make a rough estimate of the strength of 
these parameters. Clearly, the cubic and higher order invariants will influence these 
results, and could also be used to fine tune the magnetization curve. Furthermore, 
experimental evidence points toward more complicated, 8- JT] or 16-sublattice 
ordered state, which is not described by this study. However, our preliminary results on 
the 16 sublattice state seem to indicate that the shape of the magnetization curve may 
remain the same even in the case of 8 and 16 sublattice ordering in certain (possibly 
relevant) cases, where the 6r parameters of the invariants of the 7^ group are replaced 
by the appropriate coupling constant of the invariants of the larger point group. 

5. Extreme quantum case 

In this paper we have concentrated on classical spins. What changes when we take into 
account the quantum nature of the spins? When the spins are ordered, we expect some 
corrections to the theory coming from quantum fluctuations (which are in fact effectively 
taken into account by the 6's [10J). However, in spin- 1/2 quantum limit, it may happen 
that spins chose to form singlets, rather than ordering magnetically (similar physics for 
spin-1 has been studied in [4J). Even though the unit cell of the ordered state remains 
the 4 site unit cell, we need to extend our theory to take into account the inversion 
symmetry / broken by singlet bonds. We therefore need to classify bond parameters 
according of the irreps R' of the cubic group Oh = <8> T d . Having done so, we 

can write an effective energy in terms of the A R/ . For zero magnetization, the possible 
valence bond patterns are shown in figurelH when 6a 2u +26e u > 3&Ti u 1f, the valence bonds 
occupy alternate tetrahedra, otherwise they form a single valence bond in tetrahedron, 
with exchange energy distributed evenly among the tetrahedra. 

Figure 4. Possible valence bond configurations which do not break translational 
symmetry. 

If Using the standard notation for the irreps of the Oh- The irreps R of the 7^ become R' = Rg of the 
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6. Conclusions 

In this paper, we have presented a complete symmetry analysis of all possible four- 
sublattice phases of a classical pyrochlore antiferromagnet in applied magnetic field, and 
determined its magnetic phase diagram for the most general four-spin interactions. This 
theory captures the essential physics of spin-lattice coupling in spinel and pyrochlore 
materials, and provides considerable insight into the magnetization plateaux observed 
in Cr spinels. In a separate publication, we will consider the new complexities which 
arise when this symmetry analysis is generalized to the full sixteen-site cubic unit cell 
of the pyrochlore lattice [THj . 
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